Interferometric tunability of the absorption 
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We propose an interferometric setup that permits to tune the quantity of radiation absorbed 
by an object illuminated by a fixed light source. The method can be used to selectively irradiate 
portions of an object based on their transmissivities or to accurately estimate the transmissivities 
from rough absorption measurements. 
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When an object is illuminated, it will absorb radia- 
tion proportionally to its absorption coefficient: Darker 
portions of the object will absorb more light than more 
transparent ones. Is there a way around this? In this pa- 
per we analyze a setup which uses classical light sources 
(i.e. coherent beams) and permits to easily tune the 
quantity of a light absorbed by an object independently 
on its transparency, by appropriately tuning an inter- 
ferometer phase. With the same setup, high efficiency 
measurements of the absorption coefficient can be per- 
formed via a feedback mechanism. The only underlying 
assumption is that the object introduces a negligible de- 
phasing into a probe beam. Since we can employ quasi- 
monochromatic light, this assumption is met in a variety 
of systems. Moreover, in the case of objects that have a 
homogeneous phase image, the dephasing can be easily 
compensated with the interferometer phase. 

Our proposal draws inspiration from the so called 
"interaction-free-measurement" setups, where a par- 
tially transparent object can be discriminated from a to- 
tally transparent one with asymptotically negligible ra- 
diation absorption Q, 0, HI El 13 • Even though such pro- 
posals were originally based on single-photon light pulses, 
analogous results have been obtained also with classical 
light UU 

The layout of the paper follows. We start by describ- 
ing the proposed interferometric setup. We show how 
the absorption peak can be tuned and we analyze the ir- 
radiation selectivity. We then give the protocol for high 
efficiency estimation of rj. We conclude by analyzing in- 
homogeneous objects, which incorporate different trans- 
missivities. Since the process does not involve any quan- 
tum effects (such as entanglement or squeezing) one could 
also analyze it in terms of a classical theory of radiation, 
instead of the quantum formalism we use here for rigor. 

THE APPARATUS 

The proposed apparatus is a modification of the exper- 
imental setup of Ref . Q . It is obtained by concatenat- 
ing a collection of N Mach-Zehnder (MZ) interferometers 



and is depicted in Fig. ^ Initially a coherent state \a) 
enters through one interferometer port (associated with 
the annihilation operator ao), and no photons enter from 
the other port (associated with the annihilation opera- 
tor bo). As shown in Fig. after each MZ, one of the 
two emerging beams (the R beam) is focused on the ob- 
ject. Then the two beams are recombined at the input 
port of the next MZ. After N of such steps, the radiation 
leaves the apparatus at the TVth interferometer outputs 
aw and &/\r. As we will show, appropriately tuning the 
interferometers phase <j> and the number N of MZs it is 
possible to choose the value of the transmissivity rj that 
will absorb the most radiation in the object. The trans- 
missivity Tj of an object is the probability that a single 
photon will pass through it or, equivalently, the percent- 
age of the transmitted intensity of an impinging coher- 
ent beam. Note that an apparatus employing a single 
MZ which is crossed N times by the light can also be 
employed, where N can be controlled by appropriately 
tilting one of the interferometer mirrors or by using 
an acousto-optics switch. 

The input-output relations of the interferometer can be 
obtained observing that when two coherent states \a n ) 
and \(3 n ) impinge, respectively, into the input ports a n 
and b n of the n-th MZ (see Fig. |2J , the corresponding 
outputs at ports a n +l an d b n+ \ are still coherent states 
of amplitudes a n +i and /3„+i, given by 
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where (f> is the interferometer phase. The output a n is di- 
rectly fed into a port of the MZ number n + 1, while the 
output b n first passes through the object and then enters 
the other port of the same MZ. Since the object absorbs 
each photon with a probability rj without introducing any 
phase factor, its action on the input coherent state \(3 n +\) 
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FIG. 1: Proposed apparatus. It consists of N Mach-Zehnder 
(MZ) interferometers concatenated so that the output ports 
of the nth MZ is fed into the input ports of the successive one 
(for the sake of clarity the first interferometer is graphically 
enhanced). All interferometers act on the radiation with the 
same phase shift (f>. The object to be irradiated is placed out- 
side the MZs and it interacts only with the R beams. Initially 
the radiation enters from the input do. After N round trips, 
it exits through the outputs ajv and bjv. 
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FIG. 2: Mach-Zehnder interferometer constituting the nth 
element in the Mach-Zehnder sequence in the apparatus of 
FigQ The first of the two 50-50 beam splitters transforms 
the input annihilation operators a n , b n into c = (a n + b n ) / 's/2 
and d = (b n — a n )/V2 respectively. The second beam splitter 
transforms the annihilation operators c and d into a n +i = 
(c + d)/V2 and bn+i =(d- c)/V2. 



can be modeled as a beam splitter with transmissivity r\ 
that couples the input radiation to a vacuum state and 
then discards one of the two outputs. As a result the 
state \P n +i) is transformed into a coherent state of re- 
duced amplitude ^/rj/3 n +i Thus, in the presence of 
the absorber, the amplitudes a n +i and /3 n +i of the coher- 



ent states at the input of the MZ interferometer number 
n + 1 is given by 
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Iterating Eq. Q N times we can express the amplitude 
of the coherent states emerging from the whole apparatus 
as 
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Some examples of such evolution are given in Fig. 
and an analytic solution can be obtained by diagonal- 
izing S(rf) |2|. The light absorbed by the object is given 
by 



= \ao\ 



(\°ln\ 



\Pn\ 2 ) = r \a \ 2 
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i.e. the input intensity |«o| 2 minus the total output in- 
tensity I a at 1 2 + |/3jv| 2 . The quantity r is a complicated 
function of N, <f> and 77 which can be explicitly computed 
from Eq. (J5J. It measures the "effective" absorption con- 
stant of the object. 

DISCUSSION 

The possibility of changing the absorption of the illu- 
minated object from its natural value 1 — 77 to an effec- 
tive value r ~ allows one to determine the presence 
of a completely opaque object (i.e. rj = 0) with only an 
asymptotically small fraction of the input radiation being 
absorbed 0, 0, El Q] . Consider Eq. © for r] = 1 (e.g. 
completely transparent object) and r\ = (e.g. com- 
pletely opaque object). In these cases simple analytical 
solutions can be obtained yielding 



a N = a e lJV * /2 cos(A^/2) 
(3 N = ia e M ' l ' /2 sin(A^>/2) , 



for 77 = 1, and 



a N = a e tN4 > /2 cos N {<f>/2) 
Pn = 0, 



(7) 
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for 77 = 0. By choosing <fi = tt/N, from Eqs. JHJ and l|10(l . 
it is immediate to see that all radiation exits from the 
&jv-port if V = 1 and that most of the radiation (asymp- 
totically all of it for N — > 00) exits from the ajv-port if 
V = I n both cases the light absorption is minimal 
(i.e. exactly null in the first case and asymptotically null 
in the second one). Nonetheless they can be discrimi- 
nated by simply looking from which interferometer ports 
(e.g. ajv or 6jv) the light emerges. 
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The possibility of controlling the effective absorption 
r of the object by changing the interferometer parame- 
ters is evident from Fig. 0] where we plot r as a function 
of the transmissivity 77 for different values of <fi (choos- 
ing again N — n/(f>): The function r exhibits a peak for 
V = Vmax which increases from rj max — to rj max ~ 1 as 
</> decreases. This effect can be explained intuitively as 
follows. For small values of </> (i.e. high values of N) little 
radiation is leaked into the R modes at every round trip 
with the exception of the case when 77 is high. On the 
contrary, for small values of N (i.e. large values of (f>) a 
larger amount of radiation is leaked into the R modes at 
every round trip, so that the absorption peak moves to 
lower values of 77. The dependence of the absorption peak 
maximum as a function of <f> and TV is depicted in Fig. 
left. This graph also shows the values of rj max that can 
be attained in practice: it can be accurately fine-tuned 
only for r] max > 0.5 since only few discrete values of r\ m ax 
are achievable for low N, whereas high r] max ~ 1 requires 
large N, which can be difficult to achieve practically. Fi- 
nally, the selectivity of the absorption, i.e. the width 
of the effective absorption curve r as a function of 77 (see 
Fig-EJ; is n °t constant when <f> is varied: The value of the 
width-at-half-maximum is smaller for absorption curves 
peaked at r\ max ~ 0, 1 and larger for r\ max ~ 0.5. In the 
limit 4> — > the r-curve becomes a very narrow spike 
peaked just below 77 = 1. 




FIG. 3: Plot of the (rescaled) output amplitudes of the nth 
MZ interferometer |a„/ao| 2 in the a- modes at L (continu- 
ous line) and \/3 n /ao\ 2 in the b- modes at R (dashed line). 
Initially all the radiation is in mode ao, but, as the evolu- 
tions progresses, more and more radiation is transferred to 
the b-modes, until (for n = n/(j>) the radiation is entirely 
transferred. Here <j> — tv/10 so that the total transfer occurs 
for n = 10 (vertical line). Left: The object is completely 
transparent (77 = 1), so that the total energy (dotted line) is 
constant; Right the object is semi-transparent (77 = .9), so 
that the total energy decreases as the evolution progresses. 



HIGH PRECISION 77-MEASUREMENTS 

Our scheme can be easily adapted to high precision es- 
timation of the absorption coefficient, starting from low- 
quality measurements of the absorption r. [The main 
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FIG. 4: Above left: Plot of the (rescaled) apparatus output 
amplitudes jajv/ao| 2 at the output ajv (continuous line) and 
|/3jv/ao| 2 at bjv (dashed line) as a function of the transmis- 
sivity 77 of the object with (j> = 7r/60. For 77 = 1 (total trans- 
parency) all the output radiation is at fojv, whereas for 77 = 
(total absorption) all the output radiation is at ajv and a small 
amount of radiation has been absorbed. The absorbed radia- 
tion is proportional to the effective absorption constant r, de- 
fined in Eq. ©, which is depicted by the dotted lines. Above 
right: Plot of the absorption r only, in the case (f> = 7r/200. 
This plot illustrates how the setup can be used to measure 
the absorption coefficient 77 by measuring r, attaining a higher 
precision than with a direct measure (in which case the ab- 
sorption is given by the dashed line): Starting from the same 
uncertainty Ar in the measurement of r with our procedure 
we can obtain a lower uncertainty At/ in the measurement of 
77 than the one, A77, obtainable from a direct measurement. 
Below: Plot of r as a function of 77 and iV = n/(f>. Notice how 
the absorption peak shifts as a function of </>. The absorption 
peak moves to higher 77 for decreasing (f>. 



idea is revealed by the lower right graph of Fig. 0J] The 
required iterative procedure is composed by the following 
steps: i) start by roughly estimating 77 through an ab- 
sorption measurement and set the interferometer phase 
so that the r curve has a steep slope corresponding to 
such value of 77; ii) perform another absorption measure- 
ment and estimate a better value of 77; iii) again tune 
the interferometer phase, and so on. Since the absorp- 
tion curve r for the values of 77 ~ 1 can be very steep, 
a very good estimate of these 77s can be achieved even 
when the measurement of r contains a large error Ar 
(see Fig.0J). Notice that the high values of 77 ~ 1 are the 
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FIG. 5: Left: Maximum rj max (circles) and average value r\ av 
(squares) of the absorption peak r of Fig. 2] as a function 
of JV = 7r/<j!>. (The maximum and the average follow dif- 
ferent evolutions because of the asymmetry in the absorption 
curves). Increasing N (i.e. decreasing (j>), the maximum in the 
absorption peak moves to higher values of rj. The graph also 
details which are the actual values of r\ that can be achieved 
through the proposed setup as a function of N. The dotted 
line is the function [(N — 1)/JV] 4 that gives a good interpo- 
lation of the peak evolution. Right: Selectivity in the irradi- 
ation as a function of the transmissivity peak. The width of 
the peaks of the dotted line in Fig. [I] is not uniform. Here 
we plot the Root Mean Square (stars) and the width at half 
maximum (circles) of the absorption curve as a function of the 
absorption curve maximum. Notice that the RMS is almost 
constant over the whole range. 

hardest to estimate reliably without strongly irradiating 
the object, since they are associated to the region of least 
transparency. 

INHOMOGENEOUS SAMPLES 

In deriving Eq. we implicitly assumed that r\ is 
constant, i.e. that we employ a spatially homogeneous 
object with uniform absorption within the waist of the 
light beams crossing it. Instead, if it has spatially inho- 
mogeneous absorption, we can still use Eq. © to describe 
the absorption of "portions" of the incoming beam. In 
fact, in the limit in which the scale of the spatial inho- 
mogeneities is much larger than the wavelength A of the 
source light, the diffraction of the propagating beam in- 
duced by these inhomogeneities can be neglected. In this 
regime the illuminating beam can be effectively decom- 
posed in independent "sub-beams" of the wavelength A 
and interacting independently with the different portions 
of the object. A typical example is when one performs 
imaging of a macroscopic object. Controlling the interfer- 
ometer parameters we can then selectively choose which 
portion will absorb the most of the input radiation by ap- 
propriately shifting the position of the absorption peak 
of r. 

Explicitly, to selectively irradiate an inhomogeneous 
sample one thus needs to: i) estimate the transmissiv- 



ity of the various portions of the sample by conven- 
tional imaging techniques, identifying the transmissivity 
fj of the region that needs irradiation; ii) tune the phase 
4> = tt/N so that the absorption is maximized for fj. 
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